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£^ ABSTRACT 

Aims. We analyze a set of full Stokes profile observations of the flaring active region NOAA 10808. The region was 
t/~) , recorded with the Vector- Spectromagnetograph ( VSM) of the SO LIS facility. The active region produced several 

successive X-class flares between 19:00 UT and 24:00 UT on September 13, 2005 and we aim to quantify transient and 
permanent changes in the magnetic field and velocity field during one of the flares, which has been fully captured. 
| Methods. The Stokes profiles were inverted using the height-dependent inversion code LILIA to analyze magnetic field 

, vector changes at the flaring site. We report multilobed asymmetric Stokes V profiles found in the ^-sunspot umbra. 

We fit the asymmetric Stokes V profiles assuming an atmosphere consisting of two components (SIR Inversions) to 
' interpret the profile shape. The results are put in context with MDI magnetograms and reconstructed RHESSI X-ray 

images. 

Results. We obtain the magnetic field vector and find signs of restructuring of the photospheric magnetic field during 
the flare close to the polarity inversion line (PIL) at the flaring site. At two locations in the umbra we encounter strong 
fields (~ 3kG), as inferred from the Stokes / profiles which, however, exhibit a low polarization signal. During the flare 
we observe in addition asymmetric Stokes V profiles at one of these sites. The asymmetric Stokes V profiles appear co- 
spatial and co-temporal with a strong apparent polarity reversal observed in MDI- magnetograms and a chromospheric 
hard X-ray source. The two-component atmosphere fits of the asymmetric Stokes profiles result in line-of-sight velocity 
differences in the range of ^ 12km/s to 14km/s between the two components in the photosphere. Another possibility 
is that local atmospheric heating is causing the observed asymmetric Stokes V profile shape. In either case our analysis 
shows that a very localized patch of ^5" in the photospheric umbra, co-spatial with a flare footpoint, exhibits a 
00 ■ sub-resolution fine structure. 
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1. Introduction these cannot penetrate deep enough in the atmosphere, 

i n sufficient number, to directly generate the heating seen 
It is generally believed that flares are the result of ([Aboudarham fc Henouxl LL986). Fast ions have a better 
the release of stored magnetic energy during a recon- penetration range, but though these are almost certain to 
nection process acco mpanied by a change of the mag- exist in large numbers in flares, they are difficult to ob- 
netic t opology (see e.g.lCar michael 1964; St urrock fc Coppl serve. In any case the better association of both sunquakes 
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d \ 1 19661 : iKopp fc Pneumanl Il976f ). Photospheric magnetic and chro mospheric /photosphe r ic heating seems t o be with 

shear and disturbance of already emerged magnetic electrons ([Fletcher et al.l [20071 : lKosovichevll2007l ). Heating 

loops by new flux emergence are suspected to play of the mid- or upper chromosphere resulting in the emission 

an important role in t he energy loading process, see of U V radiation that then causes photospheric 'backwarm- 

e.g. jKusano et al.l (l2002l).lEngell et"ail ([201 lh or the recent ing' (jMachado et al.|[l989h is also likely. Energy may also be 

review bv lShibata fc Magaral (|2Q1 lh . carried by magnetic disturbances, whic h have the potential 

The back reaction of a large energy release on the lower r ff deep into the atmosphere (|Fletcher ^ Hnds orJ 

atmosphere has been obs erved a s photospheric 'su n quake ' ™D' but their dlssl P atl0n and heatin g 18 as Y et unknown. 

(|Kosovichev fc Zharkoval 1 19981 : iDonea fc Lindsevl 120051 ), One of the major obstacles in the study of the con- 

non-reversing change s in t h e line-of-sight magneti c field nection of photospheric magnetic footpoints and flare 

([Cameron fc Sammisl 1 19991 : ISudol fc Harvevl l2QQ5h . and events, is the lack of observations of the full magnetic 

heating of the chro mospheric temperature minimum and vector in the lower photosphere during a flare. So far 

upper photosphere ([Machado fc Rustl Il974t iMetcalf et al.l only a few full Stokes profile measureme n ts wer e recorded 

1990; IXu et al.ll2004f ). It is clear that the effects of the enor- covering a flare e.g.. iHirzberger et al.l (j2009f ) (Swedish 



mous, yet localized, release of stored magnetic energy in the Solar Telescope). Ikubo et al.l ([20071 ) (HINODE) and re- 
corona has dramatic impacts on the atmosphere beneath, centlv IWang et al.l (j2Q12h (Solar Dynamics Observatory). 
but the manner of the observed response of the deep at- Covering a larger flare site as is often required in M- and X- 
mosphere remains a puzzle. The primary agents for carry- class flares one suffers some penalty in either spatial, spec- 
ing energy in flares are thought to be fast electrons, but tral or temporal resolution. Only with sufficient spectral 
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Fig. 1. Active region NOAA 10808 on September 13, 2005. 
Left: Magnetogram obtained from SOLIS- VSM Stokes V 
profile amplitudes in the Fel 630.25 nm line. The white 
areas denote positive magnetic flux, whereas the black areas 
are regions exhibiting negative magnetic flux. The orange 
line shows the intricate polarity inversion line (PIL). Right: 
Co-spatial intensity image constructed from the Stokes / 
continuum intensities. 



sampling in the Stokes profiles can one infer sub-resolution 
structures or draw conclusions on strong gradients in 
the atmospheric parameters from the shape of the pro- 
files (jFrutiger fc Solankilll998l: ISigwarthll200ll : iMiiller et al.l 
2006; Viticchie & Sanchez Almeida 2011). For example, 
Kleint ([20121 ), revealed a complex multicomponent atmo- 
sphere around flare footpoints by analyzing high resolution 
IBIS (Interferometric B I dimensional Spectrometer) data. 

The rare full Stokes data set presented here gives us the 
opportunity to investigate the magnetic field vector changes 
associated with a flare, but due to the high spectral resolu- 
tion, also allow a more detailed analysis of the Stokes profile 
shapes. This type of detailed analysis offers the possibility 
of understanding not only the flare magnetic field evolution 
at the photospheric level, but also heating and flows in the 
deep atmosphere. 

The observed active region, NOAA 10808, showed a 
complex structure and is a so-called 5-spot, as its two oppo- 
site polarity umbrae share a common penumbra. Figure [T] 
shows the magnetogram and the continuum intensity of the 
active region. The intricate polarity inversion line (PIL) 
indicates the high complexity of the magnetic configura- 
tion typical for a J-sunspot. NOAA 10808 displayed a 
high flare activity, producing over 10 X-class flares and 
many M-class and C-class flares during its presence on 
the observable solar disk. On September 13, 2005 the ac- 
tive region was located almost at the center of the so- 
lar disk and the Geostationary Operational Environmental 
Satellites (GOES) recorded a succession of soft X-ray 
(SXR) flux peaks within 4.5 hours. Figure [2] shows the 
SXR flux registered by GOES and indicates the SOLIS 
(Synoptic Optical Long-term Investigations of the Sun) 
VSM observation time, as well as the observation time of 
RHESSI (Reuven Ramaty High Energy Solar Spectroscopic 
Imager). MDI (Michelson Doppler Imager) magnetograms 
are available for the e ntire time sectio n. The flare num- 
bering is adopted from iLiu et all ([2009) who analyzed the 
flares I - IV using a wide range of wavelength bands. 
They were particularly interested in the evolution of suc- 
cessive flares and concluded that the events were inter- 
connected. They identify the delta spot umbra as the lo- 
cation of flare III. They support their view by analyzing 
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Fig. 2. GOES X-ray flux recorde d on September 13, 2005. 
The flare numbers, adapted from lLiu et al.l (|2009f ). are in- 
dicated with capital roman numbers. The second flare is 
not discernible in the soft X-ray flux. The SOLIS- VSM ob- 
servation time from 19:36 UT to 20:57 UT is represented 
by the red solid line. MDI, recording magnetograms in Ni I 
676.77 nm and the RHESSI observing times are indicated 
with green dashed-dotted and blue dashed lines respec- 
tively. 



RHESSI hard X-ray (HXR) reconstructed images. Flare 
II cannot b e ide ntified in the GOES SXR flux. However, 
I Wang et~aH (|2007t ) observed a clear flare signature of two 
Ha kernels moving apart associated with a filament erup- 
tion. There have bee n in addition several stud i es of the indi- 
vidual flar es by e.g., [ Nagashi ma et al.l (|2007l ) , IChifor et al.l 
(2007|) and[C anou et all ([2009D as well as comparisons with 
simulation s of the large-scale emergence of a twisted flux 
rope (jArchontis fc Hoodll20Tol ). 

2. Observations 

2.1. Data reduction 

The Stokes profiles were recorded by the 
ground-based SOLI S Vector- Spectromagnetograph 
dKeller fc SOLIS Team! l2003h on September 13, 2005. 
The SOLIS-VSM has a polarization sensitivity of a few 
10 -4 per pixel in less than one second and is capable of 
obtaining the full Stokes profiles for a limb-to-limb slice 
of ~337" in width at a cadence of 5 minutes. The scans 
are performed by slit scanning with the slit positioned 
in the geocentric West-East direction and a stepwidth of 
1.125 " . The pixel size is 1.125 " and the spectral window 
encompasses the two Fe I lines at 630.15 nm and 630.25 nm 
with a resolution of 0.0027 nm per pixel. The dark currents 
are subtracted from the raw data and a flat field correction 
is employed. The polarization demodulation matrix is 
applied to obtain Stokes I, Q, U and V and fringes are 
removed. All profiles are normalized to the continuum of 



C.E.Fischer et al.: Unusual Stokes profiles 



3 



the quiet Sun intensity, which is obtained by averaging the 
continuum values of 10 quiet Sun pixel positions around 
disk center. We correlate the magnetograms thus obtained 
to counteract the spatial drift present between successive 
active region scans. Unfortunately the measurement of the 
polarized spectrum was in some parts severely affected 
by the deterioration of the polarization modulator. The 
ferro-electric liquid crystal (FLC) polarization modulator 
package is located behind the entrance slit. The deterio- 
rating effect of the FLC is the not uncommon formation 
of limited areas that do not show the expected half-wave 
retardation and/or fast axis orientation change. While 
the polarization calibration corrects for these effects, it 
cannot compensate for the loss in polarimetric efficiency. 
Therefore, the retrieved Stokes vector in the affected areas 
is still correct, but the noise is amplified dramatically. As 
the FLC is close to the entrance slit, the affected area is 
also limited in the spatial direction along the slit. The 
profiles in these regions (eastern part of the active region) 
remained noisy throughout the time sequence. Fortunately 
large parts of the umbra of the (5-spot are unaffected. 

2.2. Inversion of entire region 

We use the Levenber g-Mar quardt algorithm based code 
LILIA (jS ocas-Navarrd r2001) to perform the inversions and 
obtain the global magnetic field vector map. LILIA assumes 
local thermodynamic equilibrium and uses a height depen- 
dent model with temperature, microturbulence, macrotur- 
bulence, magnetic field strength, line-of-sight velocity, mag- 
netic inclination and magnetic azimuth as the atmospheric 
fit parameters. A (non-polarized) stray-light profile can be 
provided to the code. This profile can represent the con- 
tribution of a non-magnetic atmospheric component in the 
resolution element or/and instrumental stray-light contam- 
inating the observation. The code then determines a fill- 
ing factor of this stray-light profile which will be from now 
on referred to as the stray-light factor. We used the aver- 
aged Stokes / profile of 10 quiet Sun pixels located close to 
the active region a s the stray- ligh t prof ile. Semi-empirical 
models obtained bv lSocas-Navarrol ([20071 ) are used as initial 
model atmospheres for the umbra, penumbra and plage pro- 
files correspondingly. We a pply the non-pote ntial magnetic 
field calculation method of lGeorgoulisI ((2005) to resolve the 
180° disambiguation in the azimuth angle of the magnetic 
field vector. In this code the vertical electric current J z is 
iteratively minimized. We use J z retrieved from the disam- 
biguation of the first time scan as an initial proxy for J z 
in the following time step and carry on in this way, always 
using the previous values for the initialization, to guarantee 
consistency between the time steps. 

2.3. Asymmetric profile fits 

As LILIA uses a height-dependent model at- 
mosphere it goes be y ond Milne-Eddington 
Inversions (|Skumanich fc Litesl Il987f ) and is therefore 
theoretically able to fit asymmetries in the Stokes V 
profiles. However, it uses only one magnetic atmospheric 
component and therefore the complicated profiles caused 
by superposition of shifted profiles and opposite polar- 
ity components in a single spatial resolution e l ement 
cannot be fitted. The SIR code ([Bellot Rubiol 120031 : 
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Fig. 3. RHESSI light curves on September 13, 2005. The 
light curves were corrected for attenuator and decimation 
state changes. The rate is shown for various energy bands 
and is an average over 6 detectors. The SO LIS observation 
time is indicated with a dashed horizontal line. 



iRuiz Cobo fc del Toro Iniestal Il992[ ) is capable of using 
two magnetic components as well as an additional stray- 
light component and we employ it to invert the multilobed 
asymmetric profiles. We can not, however, achieve a 
satisfactory fit for all of the asymmetric profiles. For these 
remaining cases we therefore do not perform an inversion 
and use instead a simple approach by constructing anti- 
symmetric Stokes V profiles consisting of two gaussians 
with opposite signs. We allow the synthetic an^-symmetric 
Stokes V profiles to be shifted with respect to each other 
and to carry an optional sign reversal and vary their line 
width and amplitude. We use these parameters in a fitting 
routine employing non-linear least squares minimization to 
determine the best fit to the observed asymmetric Stokes 
V profiles. 

2.4. Context data 

We retrieve data from RHESSI and from MDI onboard 
SOHO to complement our analysis. The MDI magne- 
tograms, taken in the Nil 676.77 nm line, have a cadence 
of lmin and a pixel size of 1.98 ". The SOLIS-VSM data 
is aligned to the MDI magnetograms by cross-correlating 
the images. RHESSI recorded from 19:54 UT and the light 
curve is shown in Fig. [3l We reconstruct a low-energy X- 
ray image in the range 12keV to 18keV in the time from 
19:57:00 to 19:57:18 UT, choosing detectors 3 to 8 and 
a pixel size of 1" . A second image is obtained between 
20:04:00 and 20:04:18 in the energy range of 25keV to 
50keV with detectors 2 to 9. Both times the CLEAN- 
algorithm ([Hurford et aLl l2002) is chosen using a maximum 
of 150 itera tions. The alignment between RHESSI and MDI 
is reliable ([Fletcher et al.l [20071 ) apart from an unknown, 
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Fig. 4. Umbral region of NOAA 10808, before and after flare III. Top The TRACE whitelight intensity image on the 
left side in the upper row was recorded at 19:36 and rotated and scaled to the SOLIS-VSM data. The green contours 
are locations where the continuum intensity is below ~ 20 % of the quiet Sun continuum. These regions marked as A and 
B are also locations of low polarization signals. The second panel in the upper row is the SOLIS magnetogram which 
was obtained from the Stokes V profiles (see main text for details) recorded around 19:56 UT. The last two panels, 
magnetic field strength in kG and inclination angle 7, are the result of LILIA inversions of the SOLIS-VSM full Stokes 
profiles recorded around 19:56 UT at optical depth log(r5oonm) = — 1- The inclination angle is shown in degrees, with 
a magnetic field vector along the line of sight pointing out of the solar surface having an inclination angle of 0° . The 
black lines are contours obtained by choosing cutoff levels of 0.002 and —0.002 for the magnetogram. The black areas in 
the inversion results are non-inverted pixels. The noise produced by the deteriorating polarization modulator propagates 
into the inversion results. Bottom The first panel shows the TRACE whitelight image at 20:10 UT. The last three panels 
in the lower row are the magnetogram and the inversions results obtained from SOLIS data recorded around 20:05 UT. 
The Stokes profiles found at locations P I to P I V are discussed in section 13.21 



but small, uncertainty in the MDI roll angle. The tim- 
ing of the MDI and TRACE images was retrieved from 
the level one data fits header. The TRACE whitelight im- 
ages have a pixel size of 0.5 " and were obtained for 19:36 
UT and 20:10 UT. We used a nearest-neighbor algorithm 
to find MDI intensity images closest in time to which the 
TRACE images were then spatially aligned to. The MDI 
images were shrunk in this process to take the different 
plate scales into account. For the SOLIS-VSM images con- 
structed from the Stokes profiles we used the observation 
start time of the first slit in the active region scan. To deter- 
mine the timing of individual profiles in the SOLIS-VSM 
scans, we add the product of the slit number times the 
observation time per slit to this star t time. The RH ESSI 
instrum ent is described in d etail in iLin et ail (l2002t), the 
MDI in lScherrer et all (|l995l ) and TRACE in lHandv et all 
(|1999[ ). 

3. Results 

3.1. Magnetic vector field changes at the delta spot 

Pr eliminary inversions of this data set were presented 
by iFischer et al.l ([2009). As the deteriorating polarization 
modulator affected a wide region in the eastern part of the 
field of view, it has not been possible to deduce global mag- 
netic field vector changes for the active region as a whole. 



We focus in our analysis on flare III, an XI. 4 flare, with 
a flare site located in the umbra. The observations of this 
flare site were fortunately less affected by the instrumental 
difficulties. Figure |4] shows a close-up view of the 5-spot um- 
bra recorded at times before and after flare III. Solar north 
is in the direction of the positive y-axis and solar east is in 
the direction of the negative x-axis. The TRACE whitelight 
image shows a light bridge originating just to the south of 
the negative polarity umbra and protruding into the umbra. 
This structure stays stable throughout the entire time se- 
quence. The areas A and B are locations of low continuum 
intensities indicating strong fields. They are also marked 
in the magnetogram and exhibit in addition low polariza- 
tion signals during the entire time sequence. The profiles at 
locations P II and P IV (marked in the lower panel magne- 
togram in Fig. 3]) exhibit multilobed asymmetric Stokes V 
profiles and are, in addition to the reference profile at P I, 
examined in detail in section [3^21 P II is located at the po- 
larity inversion line (PIL) showing multilobed asymmetric 
profiles throughout the entire time sequence whereas the 
profiles at P IV, located in area A, only show asymmet- 
ric Stokes V profiles during the rising phase of the flare. 
The magnetograms were obtained from Stokes V. Instead 
of computing Stokes V at a wavelength position in the wing 
of the spectral line, we integrate the Stokes V signal over 
A. An antisymmetric spectral mask is first applied to the 
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Fig. 5. Top row TRACE whitelight im- 
age. Purple contours show areas with 
changes in the magnetic field azimuth 
larger than 5° and yellow contours are 
areas in which the longitudinal magnetic 
field switches from a positive sign to a neg- 
ative sign or the other way around. All val- 
ues are retrieved through the inversion and 
are selected at optical depth logfYsoonm) = 
— 1 . The second and third panels show the 
change in 5trans and B\ ong . White areas 
are regions of increase in Btrans/ B\ ong and 
areas of decrease are black. Rl (red con- 
tour), R2 (purple contour) and R3 (yel- 
low contour) in the last panel are loca- 
tions of simultaneous change in Btrans and 
^long greater than 400 G in both param- 
eters. Last three rows Average inclination 
angle 7 and the average i?trans i n the re_ 
gions Rl, R2 and R3. As in Fig. Hthe in- 
clination angle runs from to 180 with 
7 = 90 representing a transversal mag- 
netic field and an angle of 7 > 90 con- 
stituting a negative longitudinal magnetic 
field (vector pointing toward the Sun). The 
error bars are the standard deviation in 7 
and Btrans over a region in the umbra out- 
side the flare site. The soft X-ray GOES 
curve is overlaid in blue. 



Stokes V signal before the integration. Stokes V signals 
with the same antisymmetry properties as the mask will 
therefore appear positive whereas the integration of the 
opposite Stokes V signal will be negative. The units cor- 
respond then to the total area of the absolute Stokes V 
signal. The third and the fourth panel in Fig. 2] show the 
inferred magnetic field strength and the inclination angle 
7 at optical depth log(r5oonm) = — 1 as retrieved with the 
inversion code LILIA. The upper left areas in the inver- 
sion results are strongly affected by the noisy data and a 
'pixelation effect' is clearly visible. 

To visualize any changes in the magnetic field vector we 
subtract the longitudinal magnetic field B\ ong , the transver- 
sal field B trans an d the disambiguated azimuth angle of the 
magnetic field vector on the solar surface before and after 
flare III (19:50 UT and 20:57 UT) and display the changes 
in Fig. [5] On the left side upper panel in Fig. [5] one can 
see that the longitudinal magnetic field changes its sign 
mainly at the PIL separating the two umbrae (yellow con- 
tour on TRACE whitelight image). We use a cross correla- 
tion algorithm to align the scans taken at different times, 
after which only a slight jitter remains visible to the naked 
eye. Although a real horizontal motion of solar features 
would affect most prominently locations with strong mag- 
netic field gradients, the fact that there are several regions 



located at the PIL in the penumbra surrounding the nega- 
tive polarity umbra showing the same polarity switch iden- 
tifies an instrumental effect e.g. a remaining pointing error 
and/or a change in seeing quality the most likely cause. The 
changes in azimuth angle (purple in the same image) are 
more diffuse. There is a larger connected patch located in 
the positive polarity umbra. The azimuth angle is measured 
counterclockwise from the solar north and as the azimuth 
increases in this area, it therefore indicates a counterclock- 
wise movement of the solar surface projection of the mag- 
netic field vector. In the third upper panel we have selected 
three regions Rl, R2 and R3 at the PIL where we find 
changes greater than 400 G in both transversal and longi- 
tudinal magnetic field. As the increase in one parameter is 
accompanied by a decrease in the other, these are areas of 
real directional change of the magnetic field vector. As one 
can see in the second and third panels of the upper row in 
Fig. [5j the previously mentioned area of azimuthal angle 
change in the umbra is co-spatial with a region in which 
the magnetic field vector becomes more horizontal and is 
part of R3. We plotted the average inclination angle 7 and 
the average B trdins m the regions Rl, R2 and R3 in the last 
three rows in Fig. [5] Location Rl, at the tip of the light- 
bridge in the negative umbra, shows a decrease in 7 and, 
as 7 is above 90 , this constitutes a decrease in longitudi- 
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Fig. 6. Quiet Sun profile and Stokes pro- 
files retrieved from the positions P I to 
P IV marked in Fig. [H The columns show 
Stokes / with the two Fe I 630 nm lines and 
two telluric oxygen lines normalized to an 
averaged quiet Sun Stokes / continuum to 
demonstrate the continuum intensity dif- 
ferences and Stokes Q, U and V profiles 
normalized to the local Stokes / contin- 
uum to show the differences in scales be- 
tween the polarization signals. The wave- 
length scale is given in A A from the Fe I 
630.15 nm rest wavelength position. The 
Full- Width-Half-Maximum is given for the 
Fel 630.15 nm line. The red curve over- 
plotted on the profiles located at PI, are 
synthesized profiles fitted to the observa- 
tions by the inversion code. The magnetic 
and velocity field values at logfYsoonm = 
— 1) corresponding to this fit are B—2.1 kG 
with an inclination angle of 7=25 ° and a 
velocity of vlos=— 1.5 km/s. 



nal magnetic field. This is accompanied by an increase in 
transversal magnetic field. The trend is visible in the entire 
time series. However, there seems to be a distinct outlier 
at 19:50 UT in region Rl. This is during the rising time 
of the soft X-ray flux (purple line). At that moment the 
field appears to become abruptly more vertical as can be 
seen also in the sudden decrease of transversal magnetic 
field. Region R2 is a patch across the entire PIL separating 
the umbrae. Here 7 and the longitudinal magnetic field in- 
crease, accompanied by a decrease in transversal magnetic 
field. In region R3 there is an overall trend in the time series 
showing an increase in 7, which corresponds to a decrease 
in B\ ong as 7 < 90 , and an increase in transversal mag- 
netic field. There is a large variation between time steps 
with a zigzag pattern, leaving the impression of a sway- 
ing magnetic field vector. Region R3 is, however, located in 
the positive polarity umbra and is in the region affected by 
the instrumental difficulties mentioned in earlier sections. 
These profiles were more noisy, difficult to invert and the 
results have to be taken with caution. 



3.2. SOLIS-VSM Stokes profiles 



In Fig. [6] we plot examples of different types of Stokes pro- 
files we encounter in the sunspot. Their locations (PI to 
P IV) are indicated in the magnetogram in Fig. |4j In addi- 
tion a quiet Sun profile (location not specified in Fig. 2J is 
shown for comparison. 



3.2.1. General profile in the sunspot at P I: 

The Stokes profiles located at P I, at the umbral edge, show 
broadened Stokes / profiles. The FWHM is almost two 
times larger than in the quiet Sun (see Fig. [6]) obviously due 
to the Zeeman splitting in the strong umbral field. There 
is no significant Q or U signal, and one observes an anti- 
symmetric Stokes V profile as would emerge from a one- 
component magnetized atmosphere with the magnetic field 
vector along the line of sight and no large gradients in the 
magnetic and velocity field. We are able to easily fit such 
profiles with the inversion code, and the synthetic profiles 
are over-plotted with red lines. 

3.2.2. Asymmetric Stokes V profile at polarity inversion line 
at PIT 

The Stokes profiles for both Fe I lines located at P II at the 
PIL, display a shallow intensity profile with a strong line 
broadening in Stokes /, a linear polarization signal in Stokes 
Q and a weak, very asymmetric but not antz- symmetric cir- 
cular polarization in Stokes V. We find several such Stokes 
V profiles located at the PIL. They maintain their asym- 
metric Stokes V profiles throughout the entire time se- 
ries. T hey can be readily explain ed as so-called crossover 
points. ISkumanich fc Litesl (|l99ll ) observed that asymmet- 
ric Stokes V profiles found at the polarity inversion lines 
of J-spots are reproducible by combining 'normal' anti- 
symmetric opposite polarity Stokes V profiles with small 
velocity differences (AA in wavelength) from either side of 
the PIL. A combination of observed profiles taken from ei- 



C.E.Fischer et al.: Unusual Stokes profiles 



7 



a) Temporal evolution of 
Stokes V at P IV 



b) SIR-Inversion of Stokes V at P IV (19:58 UT) 
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c) Gaussian fit of Stokes V also in patch A (19:58 UT) 
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Fig. 7. (a) Temporal evolution of Stokes V profile of pixel P IV located in patch A. (b) 2-magnetic component SIR 
Inversion of Stokes profiles at P IV at 19:58 UT. The wavelength scale is given in A A from the Fel 630.15 nm line rest 
wavelength position. The upper panel shows the observed Stokes / profile and the Stokes / profile retrieved from the 
inversion (red dashed line). In the lower panels on the right side we plot the observed Stokes V profile (black solid line) 
and the Stokes V profile from the best-fit inversion (red dashed line). The rest wavelengths for both Fel 630.15 nm and 
Fel 630.25 nm are indicated by vertical lines. On the left side the Stokes V profiles synthesized from the two individual 
magnetic atmospheres are shown in green and blue (without taking their respective filling factors into account). The 
atmospheric parameters from this inversion can be found in Table [TJ For this fit the inversion resulted in a stray-light 
percentage of ~ 14.2 %. The line-of sight velocity difference between the two components is ~ 12 km/s. (c) Simple gaussian 
fit to the observed Stokes V profile of a pixel close to P IV and also located in patch A. The slit scan containing this 
pixel was taken approximately Is earlier than the profile at PIV. The rest wavelengths for both Fel 630.15 nm and 
Fel 630.25 nm are indicated by vertical lines. On the left side we show two anti-symmetric Stokes V profiles which 
were constructed with gaussians and shifted with respect to each other. The velocities in the two components differ 
by ~ 13.3 km/s. On the right the combination of the two synthetic an^z-symmetric Stokes V profiles is shown as a red 
dashed line. Lower row: On the left side we show a synthetic negative polarity Stokes V profile (green line) and a reversed 
polarity Stoke V profile (blue). On the right side we show again the combination of the two synthetic an^-symmetric 
Stokes V profiles as a red dashed line. The velocities in the two components differ by ~ 0.48 km/s in this case. 



ther side of the PIL is qualitatively capable of reproducing the shape of the asymmetric profile at P II. In addition the 
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Stokes Q and U profiles stay the same in magnitude and 
shape on either side of the PIL. 

3.2.3. Profiles in patch A and B in the negative and positive 
polarity umbra (PHI and PIV): 

The profiles at P III, and P IV, representing the pixels lo- 
cated in patch B and A, respectively, show a low continuum 
intensity and broadened Stokes / profiles, as expected in ar- 
eas of strong magnetic field. However, there is no linear po- 
larization visible, and the Stokes V profile is very weak (less 
than 30% than in its surrounding). Possible explanations 
for such a low polarization signal are strong magnetic fields 
canceling each other within the resolution element or /and 
a high (possibly polarized) stray-light factor or /and local 
heating causing a change in the line profile. The continuum 
intensity is at 20 % of the quiet Sun intensity at these lo- 
cations. Previous studies of the umbral atmosphere at very 
strong magnetic field locations wo uld suggest an intensity 
of rsj 10 % in our wavelength range (|Maltbvl [T970). On close 
inspection one finds that the splitting is more discernible 
in the Stokes / profile of the Fe I 630.25 nm line at P I than 
in P IV or PHI and one can differentiate between three 
distinct components in the Stokes / profile. These two ob- 
servations hint at a high stray-light factor in this region. 

The low polarization signature remains the same for 
the entire observing time in patch B. In patch A however, 
the polarization signal becomes even more reduced and we 
observe multilobed profiles in one of the Fel lines during 
flaring activity. The temporal evolution of the Stokes V 
profiles located at PIV (patch A) is shown in Fig. [Tfa). 
Between 19:47 and 20:07 UT, during the rising phase in 
the SXR flux of flare III one can observe asymmetric mul- 
tilobed profiles only in the Fel 630.15 nm line at PIV. 
The fact that only the Fel 630.15 nm line is affected is 
due to the line formation region difference between the 
two lin es that is estimated (f or non-magnetized regions) at 
64 km (jFaurobert et alJl20Q9h with the Fe I 630 . 1 5 nm being 
formed in higher layers in the atmosphere. This means that 
the atmospheric parameters responsible for the observed 
asymmetry change their properties on a very short height 
range along the line of sight in the line formation region. 
If the low polarization signal and the shortlived multilobe 
profiles would be due to a real opposite polarity emerging 
we would expect the line with a formation region lower in 
the atmosphere (Fel 630.25 nm) to be more affected and as 
this is not the case we exclude this possibility. 

We perform SIR inversions on the multilobed profiles 
taking only the spectral window of the Fel 630.15 nm 
line into account. We allow for two magnetized atmo- 
spheres within the resolution element and an unpolar- 
ized stray-light component. The best-fit synthesized pro- 
file for the Stokes V and / profile at location P IV are 
shown in Fig. [7(b) and the corresponding retrieved at- 
mospheric parameters are listed in Tab. [TJ The Stokes Q 
and Stokes U signals were below the noise. The inversion 
results in a stray-light factor of 14.2%. The two magne- 
tized atmospheres exhibit a line-of-sight velocity difference 
of — 12 km/s with an upflow of 4.64km/s and a downflow 
of 7.35 km/s. We were not able to invert some of the pix- 
els in patch A and for this cases we used a simple fitting 
approach employing constructed gaussians as described in 
section 12.31 



Table 1 . Atmospheric parameters for the two magnetic at- 
mospheres from the SIR inversion: line-of-sight velocity v, 
magnetic field B, temperature T, inclination angle 7 and 
filling factor /. The Stokes V fit is shown in Fig. [71 (b) 
with the Stokes V profile in green being the profile syn- 
thesized from umbra model 1 and the Stokes V profile in 
blue emerging from umbra model 2. Negative line-of-sight 
velocity values correspond to an upflow. 



Model 


v [kms 1 ] 


B[G] 


T[K] 


7[°] 


f 


Umbra 1 


7.35 


1050 


4046 


166 


0.2 


Umbra 2 


-4.64 


1392 


3124 


178 


0.8 



In Fig. Wis) we snow such a fit for a pixel also located in 
Patch A and close to P IV. Depending on the start model 
(same polarity or opposite polarity atmospheres) we ar- 
rive at different fit solutions and obtain two possible atmo- 
spheric scenarios that could cause the asymmetric Stokes V 
profile. The profile can be reproduced by two same polarity 
synthetic Stokes V profiles with a high velocity difference. 
The best fit achieved in this way produces a line-of-sight ve- 
locity difference of ~ 13.3 km/s. Calibrating the wavelength 
scale with the telluric line position, the blue-shifted line has 
a velocity of ~ 8.2 km/s and the red-shifted exhibits a veloc- 
ity of ~ 5.1 km/s. The second possibility is a locally heated 
atmosphere causing the complex emer ging Stokes V profile 
as has been found previously by e.g. ISocas- Navarro et al.l 
(|2000b[ ) in their umbral flash studies using the Call in- 
frared triplet. They observed shortlived asymmetric Stokes 

V profiles in the sunspot umbra and were able to show that 
these are caused by the combination of two unresolved at- 
mospheric components, one resulting in a normal Stokes 

V profile and the second producing a blueshifted reversed 
Stokes V signal due to an emission process taking place. In 
the second row in Fig. [7(c) we demonstrate this possibility 
with our simple gaussian model. We fit the profile using 
a Stokes V signal as it would arise due to absorption and 
a reversed polarity signal caused by an emission feature. 
The reversed Stokes V profile exhibits a narrower profile 
which could be due to the emission feature being present 
mainly in the core of the spectral line. The line-of-sight ve- 
locity difference between the two an^-symmetric profiles is 
-0.48 km/s. 

As the multilobed profiles appear in a strong magnetic 
field location with a high stray-light percentage in which 
the polarization signal is already reduced in amplitude, 
the effect of the spectral line going into emission is visi- 
ble in Stokes V whereas the line core reversal can be com- 
pletely masked in Stokes I as was demonstrated in IVitasI 
(|201lt). They synthesized profiles from MURaM snapshots 
(|Vogler fe Sc hiissler 2003, see code description) and showed 
that the expected line reversal in the Stokes / profile emerg- 
ing from an atmosphere with a temperature increase with 
height was masked by adding a small stray-light contami- 
nation whereas the Stokes V profiles maintained a complex 
shape. 

In both cases, large velocity difference and local heating 
scenario, a good fit is achieved using two an^-symmetric 
Stokes V profiles. This indicates a multicomponent mag- 
netic atmosphere and sub-resolution fine structure. 
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Fig. 8. MDI magnetograms during the flares I, III and V. The images are all clipped to the same maximum and minimum 
of ±500 G. The observation times are indicated by the vertical lines in Fig. [9] corresponding to the letters a to i. Upper 
row [djbjc]: The first panel shows the sunspot before flare L Panel two was taken during the rising phase of flare I and 
panel c after the flare SXR peak. Middle row [d,ej]: Magnetograms before flare III, during the rising phase of the flare 
and after. Lower row: [g,h,i] Magnetograms snowing the sunspot this time before flare V, during the rising phase of the 
flare and again after. 



3.3. Sign-reversal in MDI magnetograms and R HESS I 
reconstruction 

We select MDI magnetograms before, during and after the 
GOES SXR flux peaks of flares I, III andV. The magne- 
tograms are shown in Fig. [8] and are labeled a to i. The 
times corresponding to these letters can be found in Fig. [9] 
During the soft X-ray flux rising times of flare I (panel b in 
Fig. [8j) one can observe a patch (co-spatial with patch A 
from Fig. [4]) exhibiting an apparent reversed polarity in the 
negative umbra, before being reduced in amplitude shortly 
afterward (panel c). During flare III the reversed polarity, 
located at patch A, is more prominent during the SXR flux 
rising phase (panel e). The last flare shows an elongated 
structure exhibiting strong sign reversal and spanning over 
the PIL connecting patches A and B. 

MDI uses the difference of dopplergrams observed in 
the Stokes I — V and I + V profiles of the Nil 676.77 nm 
line to quantify the magnetic field instead of Stokes V am- 
plitude maxima as m easured by a classical magnetograph 
([Scherrer et al.lfl995f ). MD/-saturation, showing up as grey 
patches in the magnetogram, is caused by the low intensities 
measured in the strong umbra rendering the on-board pro - 
cessing algorithm ineffective as shown by iLiu et aD (|2007l ) . 
A polarity reversal, observed as a parasitic opposite polarity 
in an otherwise unipolar field however hints at a real phys- 
ical meaning. If actual flux emergence is excluded, other 



possibilities are supersonic speeds, pushing the spectral line 
partly out of the analysis window or line rev ersal due to the 
line going into emission (|Qiu fc Garyil2003l ). 

To investigate the timing relationship between the re- 
versed polarity, multilobed Stokes V profiles observed with 
SOLIS and the flare excitation we plot the total intensity 
in the pixels exhibiting the reversed polarity and the num- 
ber of multilobed Stokes V profiles as a function of time, 
and compare with the GOES SXR flux, and its time deriva- 
tive. The latter can usually be taken as a good proxy for 
the HXR time profil e (via the 'Neupert Effect', see e.g., 
iDennis fc Zarrofl993f ) which was not continuously available 
throughout the event. As shown in Fig.[9]the timing of this 
reversed polarity and the unusual Stokes V profiles coin- 
cides very well with the derivative of the soft X-ray flux. 
The appearance of the apparent reversed polarity and the 
multilobed Stokes V profiles is therefore very likely due to 
impulsive-phase, non-thermal processes which must, there- 
fore, be affecting the photospheric layers. 

This is supported by the RHESSI reconstructed images 
shown in Fig.[l0l The data was taken during the rising SXR 
flux phase (see section [2] for processing details). In our re- 
const ructions (in accordance with the findings of ILiu et al.l 
2009) we can identify a low-energy X-ray source over the 
negative polarity umbra and stretching on one side over the 
PIL into the positive umbra and at least two non-thermal 
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Fig. 9. Soft X-ray flux recorded by GOES between 19:00 UT and 23:40 UT (solid line). The dashed green line corresponds 
to its derivative. The dashed red line corresponds to the sum of parasitic positive magnetogram values found in patch 
A in the negative polarity umbra of the MDI magnetograms. The function was multiplied by a factor for the purpose 
of shifting the values into the plotting range. The blue dashed-dotted line is the number of multilobed Stokes V profiles 
observed with SOLIS and found in patch A again scaled to the plotting window. The vertical lines at the top of the 
Figure signify the record time of the magnetograms (a-i) displayed in Fig. [8j 



HXR sources placed at the location of patches A and B. 
In a classic flare picture the SXR is a loop top source of 
flare III with the hard X-ray sources located at the chro- 
mospheric footpoints of the interacting loops. These foot- 
points are areas with high-energy electrons impacting the 
lower atmosphere. 



4. Discussion and Conclusion 

We applied a height-dependent inversion to the full Stokes 
profiles data of flare III and find a change of magnetic field 
vector direction along the polarity inversion line. We iden- 
tify regions Rl and R3 which show an increase in transver- 
sal magnetic field accompanied by a decrease in longitu- 
dinal magnetic field and a region R2 in between showing 
the opposite movement (see Fig. [5]). These regions are lo- 
cated close to the PIL separating t he umbrae . With their 
extensive mult i wavelength analysis iLiu et all ((2009) have 
speculated on the flaring loop configuration (see their Fig. 
7 panel d) of flare III. Two of the flare kernels (their kl 
and k2, corresponding to our patch A and B) are in the 
positive and in the negative polarity umbra respectively 
and are connected to kernels outside the ^-spot. According 
to their speculation the flare is a consequence of reconnec- 
tion taking place between these loop systems just above 
the PIL separating the two umbrae. After the flare kl 
and k2 are connected by a low lying loop. This scenario 
would fit well with our observations, as the transverse mag- 
netic field increases on the opposite sides of the PIL (re- 
gions Rl and R3 in Fig. [5j) indicating the development 
of a lower lying loop. The increase in longitudinal mag- 
netic field at location R2 could be attributed to the re- 
moval of opposing longitudinal fields close to the X-point 
above the PIL. This is a favorab l e two -loop configura- 
tion in the flare model of iMelrosel ([19971 ) where a longer 
loop and a low-lying shorter loop are formed as a conse- 
quence of reconnect ion between two like-polarity, collinear, 
current-carrying loops crossing each other above the pho- 
tosphere. Permanent changes of the magnetic field have 
been quantified in many flares. ISudol fc Harvevl ([2005D used 



GONG magnetograms to identify magnetic field changes 
and found that the changes were often abrupt (less than 
10 m inutes) and perma nent. Also more recent observations 
from lWang et all ([2012D taken with the SDO/HMI instru- 
ment show an enhancement in horizontal field by 30 % tak- 
ing place within ~ 30 minutes in an area close to the PIL 
during an X2.2 flare. Our findings show that the changes 
occur during almost the entire time sequence (90 minutes) 
and are therefore not sudden. In addition flare III occurs 
after two previous flares, and we do not have enough cover- 
age before and after the flares to compare the field changes 
in the same area when not affected by a flare. 

We find two locations of strong magnetic field A and B, 
as inferred from the corresponding Stokes I profiles, that 
show a low polarization signal (less than 30 % than in their 
surrounding). Profiles located in patch A, such as P IV, also 
show multilobed asymmetric Stokes V profiles during flare 
activity co-spatial and co-temporal with a reversed polar- 
ity signal in MDI magnetograms and a hard X-ray source 
observed with RHESSI. The asymmetry seen in Stokes V 
lasts around 15 to 20 minutes and occurs during the impul- 
sive phase of flare III and only in the Fe I 630.15 nm line. As 
the Fe I 630.25 nm line is not affected and is formed around 
64 km deeper than the Fel 630.15 nm line, this indicates 
a very rapid variation of the atmospheric properties with 
depth. 

The temporary reversed magnetic polarity patch, if in- 
deed due to an actual magnetic feature rather than the 
disturbance of the spectral lines caused by flare heat- 
ing , is reminiscent o f the magnetic configuration described 
by iFalconer et al.l ((2000). In this, an isolated patch of 
one 'included' polarity surrounded by the opposite po- 
larity gives rise to a coronal null and a 'spine-fan' mag- 
netic t opology - whi c h is also c entral to flare s stud - 
ied by lAulanier et~aTI (|2QQ0h and iFletcher et"aD (|200lh . 
However, unlike in these models the reversed polarity patch 
in our flare is not stable, and does not persist after the flare, 
but is a transient feature. Also, due to the fact that the mul- 
tilobe profile is observed only in the Fel 630.15 nm line it 
is also reasonable to exclude flux emerging at these times. 
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Fig. 10. RHESSI reconstructions during 
flare III. The times refer to the MDI mag- 
netogram. For the RHESSI reconstruc- 
tion times refer to section [2~4l Upper row: 
RHESSI low-energy X-ray reconstruction 
(12-18 keV) contoured at a level of 40 % 
and overlaid onto a MDI magnetogram 
with the same cutoff level of ±500 G as 
m Fig. The low -energy X-ray emitting 
region is located above the PIL at the 5- 
spot. The blue box is the FOV chosen for 
the RHESSI reconstruction. Lower row: 
RHESSI hard X-ray reconstruction (25- 
50 keV) contoured at a level of 40 % and 
overlaid onto a MDI magnetogram. The 
blue box is again the FOV chosen for the 
RHESSI reconstruction. The regions de- 
noted A and B are the same as in Fig. 2] 



Resulting from our fits we find at least two possible 
atmospheric scenarios that are capable of reproducing these 
profiles. One is a two component atmosphere exhibiting a 
velocity difference between the components in the range 
of 12km/s to 14km/s or alternatively there is a polarity 
reversed contribution to the emerging Stokes V profile due 
to a temperature increase from the local back- warming from 
the chromosphere. As this effect can be masked in Stokes 
/ due to a small filling factor of the reversed signal and a 
high stray-light factor at this location we were not able to 
distinguish between these two cases. 

Observation of high velocities deduced from asym- 
met ric Stokes V profiles in a su nspot were also found 
bylMartinez Pillet et al.l (|1994| ) and lMeunier fc Kosovichevl 

(2003). In both cases th ey were observed c lose to 

the PIL in a (5-sunspot. iMartinez Pillet et al.l ([1994D 
Stokes V fit showed velocities of 14km/s. However 
th ese were long- lasting flow s and in the observations 
of lMartinez Pillet et al.l (|l994f ) seemingly unconnected to a 
flare. A transient flow of 2km/s at flare footpoints was re- 



ported recently by Ma rtinez O liver os et al.l (|201ll ) deduced 
from line shifts in the intensity profiles observed with the 
HMI onboard SDO during a whitelight flare. A downflow 
of 8 km/s was found in the N a Di line (lower chromosphere) 
bv lDonea fc Lindsevl ((2005), in association with the launch 
of a flare seismic wave, and transient dow nflows have also 
been reported earlier by iKosovichevI (|2006l ). Indeed, a pop- 
ular model for launching the flare seismic waves is strong 
local heating in the chromosphere leading t o downward- 
propa gating hydrodynamic shocks as found bv lFisher et al.l 
(1985), which deliver energy and momentum to the photo- 
sphere. 

Apparent flux reversals in magn etograms during flares 
have been reported e arlier by, e.g ., [Patterson! (|l984l ) and 
lOiu fc Garvl (|2003f ). lOiu fc Garvf (|2QQ3h also found co- 
spatial HXR sources at the reversed polarity sites. They 
simulated the MDI signal in Nil 676.77 nm and suggest 
that the signal is caused by the line center or even the whole 
line going into emission. They find an interesting connection 
between locations of strong magnetic field and reversed po- 
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larity sites and point out that this could be due to the fact 
that only at these locations is the energy requirement to 
produce the effect low enough. Our case supports this as we 
find the asymmetric Stokes profiles only at the site of lowest 
continuum intensity in the umbra (see Fig.Hlgreen contours 
on TRA CE whitelight image) which t herefo re favors the lo- 
cal heating scenario. iMaurya et al.l (|2Q12f ) have observed 
an apparent reversed polarity in HMI magnetograms dur- 
ing a 2.2 X-flare. They obtained the spectral line profile 
(6 spectral points) and find an intensity enhancement near 
the line core correlated temporally to the impulsive phase 
of the flare. Similar observational signatures are seen there- 
fore in several spectral lines from Nil 676.77 nm(MDI) to 
Fel 617.33 nm (HMI) and Fel 630.15 nm (SOLIS-VSM). 

Our analysis confirms the direct connection between the 
chromospheric flare response and the photospheric atmo- 
sphere as seen in the Stokes profiles. Although we cannot 
differentiate between the two physical explanations, it is 
clear that the flare impact sites show a fine structure in the 
photosphere. Although a fine structure in the umbra, as we 
observ ed, seems surprising, this has been already suggested 
by e.g. JSocas-Navarro et al . (2000a) who studied asymmet- 
ric Stokes V profiles in the infrared Call lines. They ob- 
served very similar short lived unusual Stokes V profiles 
and were only able to invert the profiles with a sub-arcsec 
multicomponent atmosphere in the sunspot umbra. With 
the sub-arcsec resolution of instruments such as HINODE 
it has been possible to actually observe this fine structu re 
directly in the chromosphere ([Socas-Navarro et al.| [2009). 

Full Stokes profile observations during a flare with high 
spatial, spectral and temporal resolution probing, in addi- 
tion to the photosphere, also the higher atmospheric lay- 
ers through multi- wavelength observations are needed to be 
able to differentiate between the possible atmospheric sce- 
narios leading to the anomalies seen in longitudinal magne- 
tograms and asymmetric Stokes V profiles. However, due 
to the current spatial resolution and the unpredictability of 
flare onset, the study of flare triggering mechanism, flare 
impact and photospheric response is still challenging. We 
have shown that fine structure is also important in areas 
of strong magnetic fields as encountered in sunspot umbrae 
and that the high resolution study of flare footpoints offers 
the possibility of understanding the coupling between the 
atmospheric layers in flare processes. 
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